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Based on the five-orbital model, we study the effect of local impurity in iron pnictides, and find 
that the interband impurity scattering is promoted by the d-orbital degree of freedom. This fact 
means that the fully-gapped sign-reversing s-wave state, which is predicted by spin fluctuation 
theories, is very fragile against impurities. In the BCS theory, only 1% impurities with intermediate 
strength induce huge pair-breaking, resulting in the large in-gap state and prominent reduction in 
Tc, contrary to the prediction based on simple orbital- less models. The present study provides a 
stringent constraint on the pairing symmetry and the electronic states in iron pnictides. 



PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Rp 

Since the discovery of high-Tc superconductors with 
FeAs layers [1], the symmetry of the superconducting 
(SC) gap has been studied very intensively. NMR studies 
had revealed that the singlet SC state is realized in iron 
pnictides [2-4]. The realized gap function is isotropic 
and band-dependent according to the penetration depth 
measurement [5] and the angle-resolved photoemission 
spectroscopy (ARPES) [6, 7]. This result is contrasting 
to high-Tc cuprates, where the nodal d-wave state is re- 
alized. The fully-gapped state is also supported by the 
rapid suppression in 1/Ti (oc T"; n ~ 6) below Tc in 
several compounds [8]. On the other hand, the relation 
1/Ti cx T"^ had been reported in different compounds 
[9, 10]. 

Figure 1 shows the Fermi surfaces (FSs) in the unfolded 
Brillouin zone [11]. Hereafter, we fix the electron number 
as 6.1 (10% electron-doped case). Then, the total density 
of states (DOS) per spin is 7V(0) = 0.66 eV"^ at the 
Fermi level. In the presence of the Coulomb interaction, 
antiferromagnetic (AF) fluctuations with Q « (tt, 0) is 
expected to emerge due to the nesting between FS1,2 
and FS3,4 [12]. Based on this fact, a fully gapped s-wave 
state with sign reversal (Ai^2 > 0, A3^4 < 0), which is 
called the s±-wave state, had been predicted theoretically 
[11, 12]. However, no conclusive evidence for the s±-wave 
state has been obtained experimentally so far. 

Historically, the study of impurity effects had offered us 
significant information in determining the pairing sym- 
metry in many superconductors. In iron pnictides, the 
SC state survives against high substitution of Fe sites 
by other element (more than 10%), like Co, Rh, Ni, Zn, 
Ru, and Ir [2, 13-20]. In Ba(Fei_:rCoa;)2As2 [19, 20], 
the AF ordered state is removed at x = 0.07, and SC 
state with Tc = 24 K appears. Tc gradually decreases 
as X increases, and it reaches zero at x — 0.17. This 
phase diagram suggests the relation — ATc ~ 2K per 1% 
Co substitution, and cannot be understood if the impu- 
rity pair-breaking gives —ATc > 20 K/% like in high-Tc 
cuprates. According to the first principle calculation. 



the impurity potential due to Co substitution for xz,yz- 

o 

orbitals is 1.52 eV, and its radius is only 1 A [21]. Also, a 
bulk superconducting state with Tc = 24 K is realized in 
Sr(Fei_a;Ir3;)2As2 for x ^ 0.25 [18]. These experimental 
facts would eliminate the possibility of nodal gap state. 




FIG. 1: Hole-pockets (FS1,2) and electron-pockets (FS3,4) in 
iron pnictides. jT^'^/j^ represents the impurity induced pair- 
hopping amplitude between (FS2, ±k') and (FSS, ±k). 

To explain these experiments, simple multiband BCS 
models with constant impurity potential ^ in the band- 
diagonal basis (a,/9 being the FSs), which we call the 
"orbital-less multiband model" , had been studied inten- 
sively. In the Born regime (I*^ <C (7rA^(0))^^), the s±- 
wave state is suppressed by interband scattering [22-24] . 
In the unitary regime, in contrast, the s±-wave state 
is very robust since the interband (intraband) scatter- 
ing is renormalized to zero (finite) if is constant and 
det{/'^} ^ [25, 26]. That is, Anderson's theorem [27] 
is recovered in the unitary regime. However, the lat- 
ter result should be re-examined since has complex 
momentum-dependence in usual multiorbital systems. 

In this letter, we present the first study of the impurity 
effect on the bulk DOS and Tc based on the five-orbital 
model given in Ref. [11]. By treating the multiorbital 
effect correctly, we reveal that the Anderson's theorem is 
seriously violated for the s±-wave state, due to the inter- 
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band scattering in Fig. 1, mainly via xz and yz orbitals. 
For this reason, only 1% impurities with a moderate or 
strong potential induce large in-gap DOS and prominent 
reduction in (-ATc - 20 K/%) in the BCS theory, 
which are comparable to those in nodal gap SC states. 
The present study suggests a reasonable possibility that 
a conventional s-wave state without sign reversal (s++- 
wave state) would be realized in dirty iron pnictides. 

Hereafter, we study the following 10 x 10 Nambu BCS 
Hamiltonian in the real d-orbital basis; {xy^ yz, zx) 



and eg (x^-y^, 



where 1^'^ '''' is the pairing potential. In the fully self- 
consistent T-matrix approximation, we solve eqs. (2)- 
(5) self-consistently. In calculating the DOS, however, 
we solve only eqs. (2)-(4) self-consistently, assuming the 
isotropic for FSl ~ 4. As we will show below, the 
reduction in Tc in the s±-wave state at nimp = 0.01 ex- 
ceeds 20 K for I ^ Wband- Then, the reduction in the 
SC gap is ~ 36 K if the relation A/Tc — 1.8 is supposed. 
Hereafter, we assume that this reduction in A had been 
included from the beginning. 



Ak 



A[ 



(1) 



where Hy. is the 5x5 hopping matrix of the five-orbital 
tight-binding model, which was introduced in Ref. [11]. 

A^'' = C^k^^k.ctC^k" singlet gap function in 

the d-orbital basis, where Ak,a is the gap function for 
FSa in the band-diagonal basis, and 11^'°' = (k;j'|k;a) 
is the 5x5 transformation unitary matrix between two 
representations. Then, the Green function is given by 



Fliiujn) -Gk(-it^«) 



(2) 



where a;„ = 7rT(2n -|- 1) is the fcrmion Matsubara fre- 
quency, Gk (i^k) is the 5x5 normal (anomalous) Green 
function, and Ek is the self-energy in the d-orbital basis. 

First of all, we calculate the impurity effect on the 
DOS in the SC state using the T-matrix approximation, 
which is reliable when the impurity concentration riimp is 
much smaller than unity. We consider the local impurity 
potential due to the substitution of Fe by other Sd ele- 
ments: In the present d-orbital basis, the impurity poten- 
tial is momentum-independent and diagonal. Then, the 
T-matrix for a single impurity, which is k-independcnt in 
the d-orbital basis, is given as 



T(ja;„) = (i -TC?ioc(it^,i)) ^T, 



(3) 



where ^ioc(«w„) = ;^ I]k ^k(«w„). Neglecting the 
crystalline- field splitting between and Cg, the impurity 
potential T is simply given as Ij i = ISj^i for 1 < j < 5, 
and Tj i = —ISj^i for 6 < j < 10. In the T-matrix ap- 
proximation, the self-energy matrix in the d-orbital basis 
is k-independent. It is given as 



(4) 



The gap function Ak in eq. (1) is given by the solution 
of the Eliashberg equation: 
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FIG. 2: (color online) Obtained DOS in the s±-wave state 
for m^p = 0.01. (a) |Ai,2,3,4| = 0.05 eV, and (b) |Ai,3,4| = 
2IA2I = 0.05 eV. 



Figure 2 shows the obtained DOS in the s±-wave 



SC state per spin for rii^ 



0.01 at T 



0, 



ing the isotropic gap in the band-diagonal basis. The 
number of k-meshes is iV = 512 x 512. The DOS 
is given by the normal Green function in cq. (2) as 
= ^ ^kIm{TrGk(c^ - tS)}. In Fig 2 (a), we put 
Ai,2 = — A3_4 — 0.05 eV. We also study the case where 
only A2 is changed to 0.025 eV in Fig 2 (b), consistently 
with ARPES measurements [6, 7]. Since these gaps are a 
few times larger than experimental values, the obtained 
impurity effect is underestimated. 

In Fig. 2, prominent in-gap state due to impurity 
pair-breaking emerges even for small value of / = 0.5, 
which means that only 1% impurity concentration in- 
duces strong pair-breaking effect in the s±-wave state. 
This result is inconsistent with the analysis in the orbital- 
less model, where pair-breaking due to unitary impurity 
is absent in the s±-wave state. Now, we explain that this 
discrepancy arises from the presence or absence of or- 
bital degree of freedom: In the band-diagonal basis, the 



3 



T-matrix in the normal state is given by the solution of 

^k,q — ^k,q + ^ -^k^p^'p^p q, (6) 
P 

where ^ = ^^^k^^q is the impurity potential in the band 
basis. Gp is the normal Green function that is band- 
diagonal. When k-dependence of C/k is small on each FS, 
becomes constant by replacing [/^'" with (J7k")keFSQ, 
which corresponds to the orbital- less model [22-26]. Us- 
ing g]'^^ = jj: J2k ^k' t^^e solution is given as 

f^'^ii- i''glJ-'i\ (7) 

which becomes band-diagonal as T^^ = ^l/5ioca ' '^"•/3 
when / ^ cxD, as far as det{/''} / 0. Thus, the pair- 
breaking due to interband scattering, which is described 
in Fig. 1, is absent for / = cx3 in the orbital- less model. 

However, the above discussion does not hold in iron 
pnictides since C/k strongly depends on k: For example, 
IC/^'^'^P - cos2 6lk and |C/f '^P - sin^^k for FSl (where 
6'k = tSLn~^ (ky / kx) — 7r/4), and xz,yz orbitals also con- 
stitute a large part of FS3,4 [11]. In fact, we have verified 
numerically that ^ has large offdiagonal elements [28] . 
For this reason, the impurity pair-breaking for s±-wave 
state is as large as that in the nodal gap state. 

As recognized in Fig. 2, the impurity pair-breaking ef- 
fect strongly depends on the sign of / except for I = oo: 
The impurity effect is stronger for / > 0, and it is 
the most prominent at / ~ -1-1. 5eV. This result origi- 
nates from the fact that the quasiparticle damping 7q, = 
nimpImT^ ^(— lO), which works as the depairing effect 
when the Anderson's theorem is not satisfied, depends 
on the sign of / in the presence of strong particle-hole 
asymmetry. [29, 30]. For the same reason, we find that 
the residual resistivity in the T-matrix approximation, 
which is given by c/pimp = (2e^//i) Tr{9_ff^/9/cj, • 
Gi,{+iO)-dH^/dkx-Gk{-iO)} for A (h/e^ ^ 6.45 m, 

o 

and c ~ 6 A is the inter-layer spacing) [29], is larger for 
positive /: Pimp(?iimp = 0.01) = 22, 14, and 10 [^ficm] 
at / = 1, 4.5, and oo [eV], respectively. For / < 0, 
Pimp('T'imp = 0.01) = 3 and 7 [/ificm] at / = — 1 and —4.5 
[eV], respectively. In Ref. [31], pimp for 1% Co doping in 
polycrystalline Ndllll exceeds 120 /iilcm, which corre- 
sponds to 30 fiflcm in single crystal, pimp in Co-doped 
Lallll is still larger. Near AF quantum-critical-point. 
Pimp can exceed the value given by the T-matrix approx- 
imation since the many-body correlations are enhanced 
around the impurity, as discussed in Ref. [29] . In under- 
doped high-Tc cuprates, for example, pi,np due to 1% 
Zn impurity reaches 100 pficm, which is 10 times of the 
residual resistivity in over-doped cuprates [29]. 

We also study the impurity effect on the s_|_+-wave 
state with Ai,3,4 = 0.05 eV and A2 = 0.025 eV in Fig. 
3. In this case, pair-breaking due to interband scattering 
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FIG. 3: (color online) Obtained DOS in the s++-wave state 
for 7 = 00. We put Ai,3,4 = 2A2 = 0.05 eV. 

is very small since the sign of A^ on all the FSs are the 
same. As nimp increases, the structure of DOS is grad- 
ually smoothened. The absence of zero-energy state is 
consistent with many NMR measurements [2-4, 9, 10]. 

Next, we calculate the impurity effect on Tc. To 
demonstrate the qualitative difference between s±-wave 
and s+-|_-wave states, we introduce the following BCS- 
type pairing interaction in the band basis [32] : 

^k,qi'^^ri, iuJm) = 5q,/3 ' 0k0q ' £,{^n)^i^rn) (8) 

where 0k = 1 for s+4.-wave, and 0k — sgnjcos /c^; cos fcy } 
for s±-wave. We also put S,{ujn) — ujjj/{uj'^^ +ujj-i), where 
UJD is the cutoff energy. V^'^ '''' in eq. (5) is given by 

Ea,/3^"''' • C^k"(C^k '")* • K^rK'^'- Then, is ob- 
tained by solving the eqs. (2)-(5) fully self-consistently. 
Here, we do not consider the mass-enhancement since 
the energy derivative of the self-energy Sj^(ia;„) = 
^I]q,m ^k,q(iw„,ia;m)Gq(za;m) vanishes for w„ 
since eq. (8) is a separate function of ujn and LUm [32] . 

To consider the superconductivity due to interband 
interaction between FSl, 2 and FS3,4, we put ga,i3{= 
gf3,a) = -gi for {a,f3) = (1,3), (1,4), ga,0 = -g2 for 
{a, (3) = (2, 3), (2, 4), and g^^p = for others. When 
ujd = 0.03, the obtained Tc's are shown in Fig. 4 for 
(a) 51 = 52 = 2 eV (Tco = 46 K at ^imp = 0) and (b) 
91 — = 3 eV (Tco = 40 K). In the present choice 
of 0k) TcO is equal for both s±- and s+_|_-wave states. 
In case (a), the obtained |A„| in the band diagonal ba- 
sis is almost isotropic and the same for a = 1 4: In 
contrast, |Ai| : IA2I : 1A3(4)| « 3 : 1 : 1.5 in case (b). 
In the ECS theory, -ATc = Tco - Tc is independent of 
Tco, and a qualitative relation —ATc "fa holds when 
the Anderson's theorem is violated [25, 32]. Since 7a 
takes the largest value at 7 ~ 1 eV [30], Tc for the s±- 
wave state vanishes at n-nnp ^ 0.005 for 7 = +1 cV. In 
the s++-wave state, Tc in (b) slowly decreases with rtimp 
with downward convex, since weak pair breaking occurs 
unless Aq ~ Ap. 

Until now, we have neglected the mass-enhancement, 
which is = m*/m w 2 by ARPES [6, 7]. Since the 
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FIG. 4: (color online) Obtained Tc for the s±-wave and s++- 
wave states as functions of nimp- (a) (?i,2 = 2 eV (Tco = 46 
K) and (b) 31 = 3ff2 = 3 eV (Tco = 40 K). 

depairing effect is renormalized by z, — ATc is reduced 
to -ATc • z (- 07a) [32]. When z^^ = 2, Tc for s±- 
wave state vanislies at riimp ~ 0.01, 0.02, and 0.066 for 
I ~ 1 eV, I = oo, and / = — I eV, respectively. Al- 
though -ATc ranges from 46 K/% to 7 K/%, Tc for s±- 
wave vanishes when Pimp ~ 20 /iricm for any J However, 
Sato et al. studied the impurity effects in polycrystalline 
Lallll and Ndllll, and found that Tc vanishes when 
Pimp ~ 3 mficm, which corresponds to ^ 750 pficm in 
single crystal [31]. Also, p in single-crystal Fe(Se,Te) just 
above Tc = 15 K exceeds 450 pficm in ref. [33], mainly 
due to elastic scattering. Then, the estimated mean-free- 

path Imfp is as short as apc-Fc — 2.8A [34] that is about 
1/10 of the coherence length in Bal22 [35]. (Note that 
Pimp cx is independent of z.) When Zmtp ~ ofg-Fc, 
even s++-state will be broken by localization [31]. 

In summary, we have studied the effect of Fe-site sub- 
stitution in iron pnictide superconductors. Due to the 
presence of orbital degree of freedom, the s-|--wave state 
is as fragile as nodal gap states against nonmagnetic im- 
purities. The critical residual resistivity for vanishing 
Tc 40 K for s±-wave state is only p-^^p 20 pficm. 
The corresponding mean-free-path is ^ 25aFo-Fo, which 
is longer than the experimental coherence length. 

Considering the robustness of superconductivity 
against impurities or randomness, the s_|_+ state would 
be a promising candidate for iron pnictide superconduc- 
tors. However, s±-wave state will become stable when 
(i) ]/] ^ 1 eV, or the potential radius is comparable to 
the lattice spacing and the large momentum scattering 
is suppressed. To reveal this possibility, we need more 
systematic first principle calculations for impurity poten- 
tials or measurements of pimp- Also, the s±-wave state 
can be stable when (ii) the d-orbital weight on the PS is 
completely modified by many-body effect, or (iii) strong 
coupling SC state like in heavily under-doped high-Tc 
cuprates is realized. Thus, it is important to study the 
many-body electronic states to clarify these possibilities. 
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